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ABSTRACT 



Orcus is an intermediate-size lOOOkm-scale Kuiper Belt Object (KBO) in 3:2 mean-motion resonance with Neptune, in an orbit 
very similar to that of Pluto. It has a water-ice dominated surface with solar-like visible colors. We present visible and near-infrared 
photometry and spectroscopy obtained with the Keck lOm-telescope (optical) and the Gemini 8m-telescope (near-infrared). We 
confirm the unambiguous detection of crystalline water ice as well as absorption in the Lljim region. These spectral properties are 
close to those observed for Pluto's larger satellite Charon, and for Plutino (208996) 2003 AZ84. Both in the visible and near-infrared 
Orcus' spectral properties appear to be homogeneous over time (and probably rotation) at the resolution available. From Hapke 
radiative transfer models involving intimate mixtures of various ices we find for the first time that ammonium (NHJ) and traces of 
ethane (C2H6), which are most probably solar irradiation products of ammonia and methane, and a mixture of methane and ammonia 
(diluted or not) are the best candidates to improve the description of the data with respect to a simple water ice mixture (Haumea 
type surface). The possible more subtle structure of the l.ljim band(s) should be investigated thoroughly in the future for Orcus and 
other intermediate size Plutinos to better understand the methane and ammonia chemistry at work, if any. We investigated the thermal 
history of Orcus with a new 3D thermal evolution model. Simulations over 4.5xl0''yrs with an input 10% porosity, bulk composition 
of 23% amorphous water ice and 77% dust (mass fraction), and cold accretion show that even with the action of long-lived radiogenic 
elements only, Orcus should have a melted core and most probably suffered a cryovolcanic event in its history which brought large 
amounts of crystalline ice to the surface. The presence of ammonia in the interior would strengthen the melting process. A surface 
layer of a few hundred meters to a few tens of kilometers of amorphous water ice survives, while most of the remaining volume 
underneath contains crystalline ice. The crystalline water ice possibly brought to the surface by a past cryovolcanic event should still 
be detectable after several billion years despite the irradiation effects, as demonstrated by recent laboratory experiments. 

Key words. Kuiper Belt - Methods: observational, data analysis, numerical — Techniques: photometric, spectroscopic — Radiative 
transfer 



1 . Introduction methane, as for Pluto. 



In the past decade, several large Kuiper Belt Objects (KBOs) 
have been discovered, unveiling the upper part of the size 
distribution of these outer Solar System minor bodies. Eris, dis- 
covered in 2003 ( Brown et al.ll2005h . is the largest object known 
to date, and its diameter in the 2500-3000km range, exceeding 
that of Pluto, motivated the discussion and the creation of a 
new class of objects in 2006: the dwarf planets. Several other 
large objects were discovered since 2003 (Sedna, Makemake, 
Haumea, etc) and focused the interest of scientists: the largest 
of them soon revealed a volatile-rich surface, dominated by 



Orc us was discovered in 2004 as one of the largest known 
KBOs (iBrown et al ] |2004 and is a peculiar object from several 
aspects. It revolves around the Sun in a Pluto-like orbit, in 3:2 
mean motion resonance with Neptune (with a semi-major axis 
of 39.16 AU, an inclination of 20.6°, an eccentricity 0.23, a 
perihelion of 30.26 AU and an aphelion at 48.06 AU). It is 
currently outbound, very close to aphelion. With i ts ~950km di- 
ameter ( Stan sberrv et aLllToOStlBrown et al."2010). it is now part 
of the intermediate-size objects, and is probably in a transition 
regime with respect to volatile surface content. Indeed, smaller 
KBOs with mostly featureless near-IR reflectance spectra 
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Table 1. Summarized prop erties of th e Orcus /Vant h binary sys- 
tem. R eferences are: m iNoll et alJ (l2008h . (2) iBrownet al] 
(l20Toh . (3):|Brown(200l. 





Parameter 


Value 


Reference 


System 


Angular separation 
Total mass 


0.26" 
6.3x102" 


(1) 
(2) 


Orbit 


Semi-major axis 

Period 

Inclination 


8980±20 km 
9.5d 

90° or 306° 


(2) 
(2) 
(2) 


Satellite 


Brightness Frac. 
Mass ratio 
Diameter 


8% 

0.5-0.03^ 
280-640km* 


(3) 
(2) 
(2) 



" Assuming an albedo ratio of resp. .5 and . 1 and equal densities 
* Assuming an albedo ratio of resp. 1. and .5 and equal densities 



seem to be totally depleted in volatiles (see iGuilbert et alJ 
120091 and references therein), while the largest KBOs show 
evidence for the presence of volatile ices such as methane and 
so metimes molecular nitrogen (hke Eris, S e dna and Makemake, 
cf. iMerlin etalJl2009t iBarucci et alJl2005t iBrown et alJl2007aL 
and references therein). A simple model of atmospheric escape 
of volatiles dedicated to th e Kuiper Belt region was computed 
bv Sch aller & Brownl ([2007b) and showed that the objects that 
are too small and too hot will most probably lose their pristine 
volatile surface inventory (such as CO, CH4 and N2) over the 
age of the Solar System, while large and cold objects could 
have retained these volatiles to the present day. Objects in an 
intermediate regime could have lost some volatiles but retained 
others, following the different loss rates at work for the various 
species. According to this model, Orcus should have lost all its 
volatile content. However, traces of methane (to be confirmed) 
could explain the near-infrared spectrum of Orcus (Barucci et al. 
I2OO8I) . This object can therefore provide key constraints on 
the current surface volatile inventory of intermediate-size KBOs. 



2.5 
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Fig. 1. A compilation of available visible spectra of Orcus, pre- 
sented at a reso l ution o f R~600. From bottom to top: this paper, 
Fornas i er et alJ (l2004 . Ide Bergh et"al. (2005), Forna sier et al.1 
720091) . The three top spectra were shifted by 0.5 in reflectance 
for clarity. 



scenario would then be an initial capture of Vanth, followed by 
large oscillations of eccentricity an inclination (Kozai cycling), 
and when the pericenter drops to a low enough value, tidal 
evolution can finally lead t h e orbi t to its current circular shape. 
As noted by iBrown et al.l (12010), a future discovery of any 
Orcus coplanar satelhte would rule out this possibility (as for 
the Pluto-Charon system). Relative color measurements by the 
same team indicate that Vanth is significantly redder than Orcus 
(another unique property among KBO binaries), a characteristic 
that is currently difficult to understand in the framework of a 
formation by a giant impact. 



So far, Orcus' rotation properties are not uniquely defined: 
0.03 to 0.18 magnitude amplitude variations are measured over 
a period ranging from ~7 to 21h (Duffard et al. 2009, and ref- 
erences therein), although I Sheppard (200 7) finds no lightcurve 
variations within the photometric uncertainties of his measure- 
ments. Such a flat lightcurve could be diagnostic of a very 
slow rotation rate, a small peak-to-peak lightcurve variation or 
a nearly pole-on geometric aspect. Also, given the size of Orcus 
and the relatively slow rotation rate, we can most probably ex- 
pect a circular shape from the relaxation of a fluid body in hy- 
drostatic equilibrium. 

To date, Orcus is known to have a single relatively large 
companion, Van th, in a clos e circular orbit (see details in Table 
[1]). According to lBrown et a l. (2010), the satellite is on a nearly 
face-on orbit (with respect to the Sun). Again, Orcus seems 
to be a unique case in an intermediate regime: larger KBOs 
have relatively small satellites in circular orbit s, while small 
KB Os have ve ry large satellites (iNoll et alJl2008l) . As described 
by iBrown et' al. (2010), two scenarios could be invoked to 
explain the creation of the Orcus binary system which we will 
summarize below. In any case, because the orbit is circular, tidal 
evolution was at work at some point of the system's history. 
First, a catastrophic collision could be at the origin of the forma- 
tion of Vanth, as hypothesized for Pluto-Charon (Canup 2005). 
Because of the large inclination of Vanth's orbit, Kozai cycling 
(iPerets & Naod 12009) could also have occurred. A second 



In this work, we review the existing spectroscopic investi- 
gations of Orcus and objects with similar spectral and physi- 
cal properties and search for additional constraints on the pres- 
ence of volatiles (and maybe their irradiation products) on an 
object that might have retained some of its original content. We 
first report additional 0.4-2.5/im photometry and spectroscopy 
obtained from the Mauna Kea large telescopes. We present 
Hapke radiative transfer modeling of the spectral data, the test- 
ing for the presence of water ice, methane and ammonia ices 
(and their irradiation products, ethane and ammonium resp.) at 
the surface and discuss the results. Finally, we use a new 3D 
thermal evolution model dedi cated to the interiors of KBOs 
dGuilbert-Lepoutre et al.ll20Tol) to describe Orcus' thermal his- 
tory over the age of the Solar System and the observable conse- 
quences of this history on the current surface. 

2. Observations and data reduction 

2.1. Visible data 

Visible spectra of Orcus were acquired in visitor mode on 
2005 April 9, at the Keck-1 telesc ope with the LRI S-RB (Low 
Resolution Imaging Spectrometer. lOke et al.l [19951) instrument 
mounted at the Cassegrain focus. Beam-splitters separate the 
light between two arms, red and blue. The field of view (FOV) 
of both arms is ~7'x6.8'. The blue side is equipped with a 
mosaic of two 2Kx4K Marconi (E2V) CCDs with a plate scale 
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Table 2. Comparison of the visible spectral slope S (in %/100 nm) as computed from linear regression using the same algorithm 
on published data and our work over the common spectral range (0.51-0.76 jum) 



Data reference 


UT date 


r(AUi;^ 


<t 


Instrument 


Computed S (%/lOOnm) 


This paper 


2005-04-09 


47.67 


0.99 


Keck + LRIS-RB 


2.01±0.54 


Fornasier et al. (2009) 


2008-02-04 


47.81 


0.48 


VLT + FORS2 


1.75±0.57 


deBerahetal. (2005) 


2004-04-11 


47.62 


1.02 


VLT + FORS2 


1.42±0.74 


fornasier et al. (2004) 


2004-02-29 


47.61 


0.46 


TNG + DOLORES 


1.80±0.67 



° Heliocentric distance from JPL Horizons ephemeris 
* Phase angle (°) from JPL Horizons ephemeris 

Table 3. Near-infrared photometry from Gemini-i-NIRI (this work) and near-infrared published colors from Orcus. 



UT date 


J ±a 


H±cr 


a:.s±o- 


J-H±cr 


H- K±cr 


Reference 


2005-02-20 


17.86±0.08 


17.73±0.05 


17.82±0.05 


0.14±0.10 


-0.09±0.07 


This paper 


2005-02-21 


17.82±0.04 


17.70±0.04 


17.81±0.04 


0.12±0.07 


-0.11 ±0.06 


This paper 


2005-02-22 


17.85±0.05 


17.74±0.05 


17.77±0.06 


0.11±0.07 


-0.03±0.08 


This paper 


Average 


17.84±0.11 


17.72+0.09 


17.80±0.08 


0.12±0.14 


-0.08±0.12 


This paper 


Aver. 2004-04-11 to 21 








0.13±0.05 


0.04±0.050 


de Berah et al. (2005) 



of 0.135"/pixel. The red arm has a Tektronik 2048Kx2048K 
detector with a scale of 0.2H"/pixel. We used the gratings 
400/3400 on the blue side and 400/8500 on the red side with 
dichroic 560. The resulting wavelength coverage is 0.37-0.83 
fiin with a resolution of ~350. We used a long slit of one 
arcsec width that we maintained along the parallactic angle 
during all the observations to minimize the loss of flux due 
to the atmospheric refraction. Arc images were obtained from 
Hg/Ne/Cd/Zn lamp frames on the blue side and Ne/Ar lamp 
frames on the red side. To avoid additional known distortion 
effects due to the telescope flexure, we acquired spectral 
flat-field and arc frames for each science telescope position. 
We corrected the spect ral images f or bias, flat-field and cosmic 
rays with the IRAF dTodvl Il986l) spectroscopy pack, which 
we also used to calibrate the wavelength. For each instrument 
arm, we extracted the spectra using IRAF/APEXTRACT 
AR4LL task and merged them into a single ID spectrum. It 
was corrected by a solar analog (HD 126053) observed close 
in airmass and time and normalized to unity around 0.55 fim. 
The resulting 2700s spectrum (which is a combination of 
three individual spectra, 900s each) is presented in Fig. [T] We 
do not believe that the slight turn-off at 0.42 yum is real. We 
inspected the three individual spectra and found that although 
the overall spectral slope remains consistent from spectrum 
to spectrum, small variations were observed near the blue 
end of the individual spectra, which may be caused by the 
poor CCD sensitivity in this wavelength region. Therefore, we 
concluded that the slight turn-down feature that appears at 0.42 
fim in the combined spectrum of Orcus is not real, but an artifact. 



We also acquired a series of V band images (through the 
broadband filter centered at 5437A and with a FWHM of 
922A) to monitor the photometry over the course of the spec- 
troscopy observations. Data reduction, photometric calibration 
an d measuremen t s were applied following a method described 
in iDelsanti et al.l (12001 ). We found an average value of V = 
19.06 ± 0.03, which coiTesponds to mv(l, 1) = y - 5 log(rA) - 
2.32 + 0.03 (with r and A the helio- and geocentric distance 
respectively in AU). Rabinowitz et al. (2007) found a value of 
Hy = 2.33 ± 0.03 from their solar phase curve study. They pub- 
lished a phase correction coefficient for Orcus of 0.114 + 0.030 



Table 4. List of near-infrared Gemini -hNIRI spectroscopy obser- 
vations. 



Gj2, UT start date Expj^ Airmass Solar An. Airm. 

2005-02-20 09:43 r20? 1.10-1.12 HIP51125 r24~ 

H 2005-03-01 10:54 2160 1.22-1.41 HIP47071 1.24 

K 2005-02-22 09:38 3600 1.09-1.13 HIP51125 1.28 

K 2005-02-22 10:52 2100 1.15-1.23 HIP51125 1.28 



« Grism. J: 1.00-1.36 yum, H: 1.43-1.96 yum, K: 1.90-2.49 fim 
* Exposure time in seconds 



mag. deg ' in V. As we observed with a phase angle of 0.99°, 
our phase-corrected absolute magnitude is therefore 2.21 + 0.04. 

2.2. Near-infrared data 

Data were acquired in queue mode at the Gemini North 
Telescope with the NIRJ (Near InfraRed Imager) instrument on 
UT 2005 February 20-22 and UT March 1, (see Table gli. The 
detector is a 1024x1024 ALADDIN InSb amy. We used the f/6 
camera, which gives a FOV of 120' x 120' and a pixel scale of 
0. 1 17 arcsec/pixel in imaging. Photometry was acquired through 
the broadband filters J, H, K short, centered at 1.25, 1.65, 2.15 
fim respectively. We used the GEMINI NIRI package to produce 
the master dark and flat-field frames. To correct for the sky 
background contribution and recombine the sub-frames into 
the final image, we used the IRAF/XDIMSUM package. We 
measured the object flux using an aperture correction method 
(see Delsanti et al. 2004, for details on the complete procedure). 
Photometric coefficients (zero-point, extinction) were obtained 
from a least-square fit for the standard star flux measurements. 
Results are presented in Table [H Our near-in frared colors are 
compatible with those of Ide Bergh et alJ (l2005h . 

We acquired the spectra through the 6-pixel wide slit and 
the f/6 J grism (1.00-1.36 fim useful range with the "blue" slit, 
resolution of -480), H grism (1.43-1.96 jum, R~520), and K 
grism (1.90-2.49 yum, R~520). Circumstances are detailed in 
Table H) We used the usual AB nodding technique to acquire 
the sub-frames of a cube of data (the telescope was dithered by 
~20" in between sub-images). Solar analogs were observed to 
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Table 5. Water ice near-infrared absorption band depths and cen- 
tral positions for Orcus. 



2.0 



2005-02-22 UT 09:38 



Data 


Central A (yum) 


Depth (%) 


H grism 


1.50 


25.5±3.6 


H grism 


1.65 


13.2±4.5 


K grism UT 09:38 


2.05 


28.2±3.6 


K grism UT 10:52 


2.05 


22.2±6.7 



correct science spectra from the telluric lines and the contribu- 
tion from the Sun's spectrum. We used the GEMINI GNIRS 
package to produce the ancillary frames (dark, flatfield, arc). 
Cubes of dithered spectra were corrected from the dark cur- 
rent, flatfield, and the distortion effects. The AB pairs were sub- 
tracted (A-B and B-A) to correct from the sky contribution and 
were re-aligned into a final composite 2D image. The ID spec- 
trum was extracted with the IRAF/APEXTRACT APALL task. 
Wavelength calibration was performed with IRAF/ONEDSPEC 
package. 



3. Results 

3.1. A featureless neutral visible spectrum 

Orcus shows a linear and featureless spectrum (within the noise 
level) over the visible range, which is compatible with previ- 
ously published data. Figure [T] shows the spectra available in 
the literature, presented at a resol ution of R~600 (data from 
this work and Fornasier et al. '2004', were degraded to match the 
R~600 resolution of the F0RS2 data). We computed a linear re- 
gression using the same algorithm for all published spectra and 
our work (on the original, non-degraded data), over the wave- 
length range that is common to all spectra (0.51 - 0.76 fim) 
in order to compare the results: values of the spectral slope in 
%/lOOnm are reported in Table |2] They are compatible at one 
sigma level and are characteristic of a quasi-neutral surface with 
solar colors. We do not believe the slight turn-off at 0.4 micron 
is real as explained in Sect. 12.11 Orcus shows very consistent 
visible spectral properties over time, which strongly points to a 
homogeneous surface, at least in the optical. 

3.2. Crystalline water ice, ammonia, and methane? 

We adjusted the visible and near-infrared parts of the spectrum 
in reflectance using the V magnitude computed in Sect. 12. II and 
near-infrared photometry from 2005 February 20, (see Table 
[3] i.e. photometry that was observed simultaneously with the J 
band spectrum). In the K band, we obtained two consecutive 
spectra on 2005 February 22, which are plotted at the resolution 
of the instrument in Fig. |2] along with the average spectrum. The 
complete reconstructed spectrum at the instrument resolutions 
from 0.40 to 2.35 fim is presented in Fig. |4] (using the average 
of K-band spectral data). 

In the near-infrared, the spectrum of Orcus shows deep ab- 
sorption bands centered at 1.5 //m, 1.65 fim, 2.05 //m (attributed 
to the presence of crystalline water ice) and smaller features 
around l.lfim. Measured water ice band depths and positions 
(using a continuum removal plus Gaussian fitting method) are 
reported in Table |5] Figure |3] shows a close-up of the data in the 
1. 4-2. 4yL/m range, the synthetic reflectance spectrum of different 
pure and hydrated ices, and some of their irradiation products 
(with 10 i-im size grains). From these, we can unambiguously 



1.5 



1 .0 




Average 



0.5 




1.9 2.0 2.1 2.2 

Wavelength (,um) 



2.3 



2.4 



Fig. 2. Orcus spectra taken with the K grism on 2005-02-22. Top 
and bottom spectra were shifted by -1-0.6 and -0.6 resp. in re- 
flectance for clarity. The solid line is a running mean to guide 
the eye (smooth box with a width of ~ 55A). 



identify crystalline water ice (with the broad 1.5 yum and 2.0 
/vm features and the 1.65 jum absorption band). Methane or 
ammonia ices ma y account for the 2. 2 jum absorption band as 
surmised also by Barucci et al.l (l2008l) . However, hydrated and 
pure ammonia ices do not share the same features around 2.2 
jum. The Orcus spectral shape beyond 2.23 //m may also suggest 
the presence of methane ice. Ethane ice is also surm ised, as it is 
an irradiation product from methane (Bar atta et alJ l2003): both 
components best described the near-infrared spectral behavior 
of Quaoar ( Schaller & Brown 2007a) . which is very similar to 
that of Orcus. We also show a spectrum of ammonium (NH^ ), as 
it is an irradiation product of ammonia and helped to reproduce 
the 2.2^m feature s in the spectrum of Pluto's satellite Charon 
(ICook et al.ll2"009l) . 



4. Radiative transfer modeling of the surface 
composition 

4. 1 . Algorithm and parameters 

To better investigate the surface properties of O rcus, w e ran a 
radiative transfer model following the theory of iHapkd (ll98lL 
1993) over the complete wavelength range available. This algo- 
rithm fits a synthetic reflectance spectrum to the data, using re- 
flectance properties (optical constants) from a user-defined input 
list of chemical components, their relative amounts, and a grain 
size. When applying this method, one has to keep in mind that 
the best-fit synthetic spectrum obtained is not a unique solution 
to describe the object's surface composition and that some pa- 
rameter values might lead to the best fit, although with no phys- 



A. Delsand et al.: A portrait of Plutino (90482) Orcus 



5 



1 .0 
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-0.5 
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1 .4 
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Wavelength (^tm) 
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Fig. 3. A close up of Orcus spectrum in H and K bands (top, grey line) and the reflectance of, from top to bottom: pure ammonia 
and, shifted by several units for clarity, pure methane ice, amorphous water ice, crystalline water ice, hydrated ammonia (H20:NH3, 
99: 1), ammonium and ethane. Synthetic ice spectra were obtained using Hapke theory (.Hapke 1 98 L . 1 993.) with the optical constants 
quoted in this paper and a particle size of 10 yum. 



ical meaning. A detailed description of this model and its limits 
are presented in lMerlin et al.l (1201 Obh . 

We chose a symmetry parameter of v = -0.4 and a back- 
scattering par ameter B - 0.5. These values are close to those 
computed bv Verbiscer & Helfensteinl(ll998l) to describe the sur- 
faces of the icy satellites of the giant planets. Finally, we com- 
puted the_geometric albedo at zero phase angle from Eq. 44 of 
[Hapke (1981). We used the Marqvardt-Levenberg algorithm to 
obtain the lower reduced value for the fit of the synthetic 
spectra to our data. Each synthetic spectrum is obtained using 
an intimate mixture of different compounds (from a user-defined 
input list). The free parameters are the grain size and the rela- 
tive amount of each chemical compound. We emphasize that the 
results presented here are not unique and only show a possible 
surface composition. 

4.2. The choice of the chemical components 

As input to the radiative transfer model, we used optical con- 
stants of several ices at low temperature (close to 40K) that are 
likely to be present on the surface of an icy body located at ~50 
AU such as: 



- Water ice jGrundv & Schmit tlll998h 

- Pure methane ice ([puirico & Schmittlll997h 

- Pure ammonia ice ( Schmitt et al.ll 19981) 

- Diluted ammonia ice (Ted Roush, personal communication) 

- Pure methanol ice (Quirico & Schmitt, personal communi- 
cation). 

We also used optical constants of a dark compound such 
as amorphous black carbon (Zubko et al. 1996), which is 
able to reproduce the low albedo of the majority of TNOs 
dStansberrv et al.l l2008h . Space weathering is known to be 



responsible for the alteration of the KBOs surfaces with time 
through the action of galactic cosmic rays ions, solar wind 
plas ma ions and energetic particles, dStrazzuUa et al.l 119911: 
Cooper et jilj |2003(). The effect on ices is documented from 
laboratory experiments (for instance [Brunetto et al.ll20()6l) : the 
fresh and bright icy mantle becomes spectrally redder and 
darker with irradiation. Without any better constraints on the 
nature of the irradiation products present on the surface of these 
objects, we decided to use tholins, which reproduce the spectral 
slope behavior best in the visible (in our case. Titan tholins were 
found to be the more appropriate; this compound is generated 
from the irradiation of mixtures of N2 and CH4, .Khare et al] 
[l984l) . 



Additional spectroscopic observations of Orcu s in the near 
infrared (H-i-K band only, from lBarucci et al.l2008h showed deep 
absorption features at 1.5 and 2.0 ;um and shallow bands at 1.65 
fim and close to 2.2 /im. The presence of crystalline water ice is 
evident in the 1.65 yum feature. The weak 2.2yum feature is much 
more difficult to interpret (owing to the detection level and mul- 
tiple candidate components). These authors tested Hapke mod- 
els of ammonia hydrate or pure methane ices (separately): none 
of these attempts could satisfactorily reproduce the shape of the 
2.2fim feature, although ammonia hydrate was the favored ex- 
planation. In this work, we will try to better reproduce the gen- 
eral shape of our spectrum on the full 0.4-2.4 yum range with 
the help of the Hapke model ing of mixtures of ices hsted above, 
assummg an albedo of 0.2 (tStansberrv et al.ll2008h . Below, we 
will not discuss in detail the various relative quantities listed in 
Table|6] because we consider Hapke model results as only qual- 
itatively indicative of the possible candidate components. The 
reasons are that: (i) synthetic models, even when accurately fit- 
ted to the data, do not provide a unique, definitive solution; (ii) 
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the SNR of the data currently available for Orcus do not allow 
this detailed interpretation; (iii) the optical constants from ma- 
terials used here where measured in conditions close to those 
that should prevail on Orcus, but some of them are missing and 
approximations were used. 

4.3. A simple water ice mixture 

Given the obvious features of water ice in its crystalline state in 
our data, our first attempt was to use a simple mixture of amor- 
phous and crystalline water ice at ~40K. Again, Titan tholins 
were used only to reproduce the spectral shape in the optical 
range (as in all models of this work): this compound is certainly 
not present as such at the surface of Orcus. Similarly, amorphous 
black carbon was only used as a mean to match the low measured 
albedo. The corresponding best-fit model is plotted in the bottom 
of Fig. 3] the quantities are tabulated in Table |6] (model 1). 

The large path lengths of amorphous water ice should be 
taken with great care. Indeed, when alternatively using a blue 
component to fit the continuum (using for instance the re- 
flectance properties of kaolinite from the ASTER database, 
http://speclib.jpl.nasa.gov/search-l/mineral), the path lengths 
and the proportions of amorphous water ice in the final model 
drastically drop. This indicates that crystalline water ice is 
most probably dominating the mixture. All previous work on 
Orcus reports the presence of water ice dFornasier et al.l 120041: 
IdeBerph etaTI 12005: Truiil lo etal]|2005h . and the most recent 
higher SN R data show the presence of crystalline water ice 
dBarucci et al.. 2008, and this work). The crystalline water ice 
fraction for this mo del is comparable to the one reported in 
iBarucci et aP (l2008h (about 15%). The depth of the absoi-ption 
feature at 1.65 yum depends on the crystalline water ice abun- 
dance and temperature of the surface. At low temperature (close 
to that of Orcus, ~40K), a deep absorption feature as observed in 
our spectrum (or in previous works) suggests a high crystalline- 
to-amorphous water ice ratio. 

This first attempt shows that a simple mixture of water ice 
approximately reproduces the general spectral behavior of Orcus 
in the near-infrared, but fails to describe the more subtle fea- 
tures (at 1 .65;um and around 2.2yum). Orcus therefore has a water 
dominated surface, but most probably hosts additional absorbing 
substances, such as volatiles, in smaller amounts. 

4.4. Adding metliane and ammonia ices 

For this second step, we tested various ices that could absorb 
around 2.2/im, contribute to the blue slope long-ward 2.3/vm 
and could possibly be present on a ~40K KBO (see also Fig. 
[3] and Sect. 13.21 ). In particular, ammonia hydrate has an absorp- 
tion feature around 2.21yum: it was detected on Uranus' satelhte 
Miranda (Bauer et al. 2002) and possibly identified on Quaoar 
( Jewitt & Luu 2004) and Orcus (previous work by Baruc ci et al.l 
^008). Methane (dominating the surface composition of large 
KBOs, such as Pluto, Eris or Makemake) is also an interesting 
candidate, with its absorption feature around 2.2/im and blue re- 
flectance in the 2.3-2.4//m range. We also included methanol in 
our tests (pure CH3OH from Quirico & Schmitt, priv. comm.) 

At this stage, best-fit results were obtained using a mixture 
of water ice (amorphous and crystalline), methane and ammo- 
nia (pure or diluted). The resulting models 2 and 3 are plotted in 
Fig.|4](see also Table|6]l. They have an equivalent reduced and 
reproduce almost equally the general shape of the near infrared 
spectrum. The agreement of the synthetic spectra to our data is 



however not totally satisfactory around 1.65jum and 2.2jum, al- 
though the slope beyond 2.3jum is now correctly mimicked. Only 
traces of methane and ammonia are sufficient to significantly im- 
prove the fit with resp ect to a simple water ice mixture. We tested 
the mixtures used bv lBarucci et al.l (l2008l) against our data, i.e. 
fi-ying separately methane, then ammonia hydrate (respectively 
mixed with water ice) but obtained a worse fit than models 2 and 
3. In particular, the 2.2yum band was poorly reproduced in both 
cases, and for the ammonia hydrate model, the synthetic spec- 
trum had a reflectance that was too high beyond 2.3/^m. 

From a theoretical point of view, ammonia and methane 
could be major volatiles in the solar nebula at the distance where 
Orcus might have formed (Le wis , 197 2). Ammonia on one hand 
is known to lower the temperature of the melting point when 
mixed with water ice. Methane on the other hand is almost insol- 
uble in aqueous liquids at low pressures (no electric dipole mo- 
ment). At the central pressure of Orcus (few kilo-bars), its solu- 
bility should nonetheless reach about 1%. Methane could there- 
fore trigger so-called explosive aqueous cryovolcan ism in the 
presence of ammonia and water dKargellll992Lll995h . This pro- 
cess has been introduced to explain the cryovolcanism observed 
on Triton dKargel & Stromlll990i) . Indeed, the bubble content of 
liquid melts would increase as they propagate toward the surface 
through cracks: the solubility of methane in these melts depend 
on the pressure, which decreases as the melts reach the surface. 
If the bubble content becomes sufficient, the melts would eventu- 
ally disintegrate into a gas-driven spray producing the explosive 
cryovolcanism. 

4.5. Searcli for daugliter species of metliane and ammonia 

Ethane is an irradiation pro duct (from UV photolysis) of solid 
ethane ("Barat ta et al.ll2003l) . and is expected to be involatile at 
~40K. Therefore, once produced by solar irradiation, it should 
remain at the surface of the object and still be observable. 

oore & HudsonI (l2003i) report that the production of C2H6 is 
more efficient when CH4 is pure. C^Hf, w as suggested at the sur- 
face of Quaoar ("Schaller & Brown 2007a') and detected on dwarf 
planet Makemake (Brown et al. 2007a). Since our data could be 
compatible with traces of methane at the surface of Orcus, we 
decided to further investig ate the contribution ofCaHf,. We used 
the optical constants from lOuirico &Schmittl(ll997h . Following 
a similar reasoning, we decided to test the presence of ammo- 
nium (NH^), which can be produced for instance by irradia- 
tion of pur e NH3 or mixtures of water and ammonia as expe- 
rienced by (i Moore et al.ll2007|) . Other species can be produced 
from the irradiation of NH3 orNH3:H20, such as hydroxylamine 
(NH2OH); unfortunately, the corresponding weak infrared fea- 
tures are not easily detectable. In this test, we used optical con- 
stants of NH^ derived by J. Cook (priv. comm.) from the studies 
by Moore et al. 

We therefore tested ammonia and methane separately along 
with their respective irradiation products. Model 4 (middle quad- 
rant of Fig. nil shows a mixture of water ice, methane, and ethane, 
while model 5 includes water ice, ammonia hydrate, and ammo- 
nium. The methane option (model 4) gives a lower reduced 
and overall better fit than model 5, and is visually and statisti- 
cally equivalent to models 2 and 3. However, ethane is finally 
found not to contribute to the final mixture (0%) when mixed 
with only CH4 and water ice, which tells us that methane is itself 
an interesting component to describe the H and K band data. The 
ammonia and ammonium mixture model is not satisfactory to 
improve the overall fit of Orcus spectrum (with respect to other 
models), but the structure of the 2.2//m band is worth of interest. 
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Fig. 4. Visible and near-infrared spectrum of Orcus (grey line) with synthetic spectra superimposed (solid black line). Spectra are 
shifted in albedo for clarity from -0.1 to about +0.5 (except model 2). Quantitative description of the models and reduced are 
presented in Table |6] The bottom spectrum shows that a water ice-only model (Haumea-like surface) fails at describing all of Orcus 
features, especially in the K band, although water ice is the likely dominant component. The two models above are the best-fit results 
from our next investigation step involving volatiles: they include a mixture of methane and ammonia (pure of diluted) in addition to 
water ice (in its amorphous and crystalline state). The middle part of the plot illustrates the models including ammonia or methane 
and their respective irradiation products (ammonium and ethane). The top part shows our overall best-fit result; a mixture of water 
ice, ammonium (NH|), and ethane (CaHe). 



4.6. Best-fit results and conclusions 

For this last step we aimed at better describing the data overall, in 
particular in the 2.2/im region, using results (and possible candi- 
dates) derived from the previous investigation steps. Our best-fit 
result overall is displayed as model 6 (top of Fig.|4|i and includes 
a mixture of water ice mostly in its crystalline state (~20%), am- 
monium (~13%) and traces of ethane (~7%). The overall shape 
of the Orcus spectrum is correctly reproduced, in particular in 
the K band with both the 2.05jum large water ice band and the 
2.2//m region more correctly fitted. The presence of ammonium 



also better reproduces the visible slope. The 1.65jum feature was 
never completely correctly reproduced, which could be related 
to the SNR of the data. Our conclusion is that the irradiation 
products NH^ and C2H6 greatly helps to describe the data in the 
1.4-2.4//m region, although both parent species, as a mixture 
of ammonia (pure or diluted) and methane cannot be firmly ex- 
cluded (model 2 and 3 of equivalent statistical level). We believe 
that the four species are probably present as traces at the surface 
of Orcus. Therefore, further spectroscopic studies of Orcus with 
significanlty higher SNR (especially in the 2.2/im region) would 
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Table 6. Quantitative description of the chemical composition models with the relative percentage of each component. Grain size 
in micron is indicated in parenthesis. Models are sorted on Fig. |4]from bottom (model 1) to top (model 6). 



Components 


Model 1 


Model 2 


Model 3 


Model 4 


Model 5 


Model 6 


Reference 


Amorphous H2O 


13% (500) 


5% (600) 


7% (400) 


9% (370) 


5% (300) 


1% (200) 


Grundv&Schmitt(1998) 


Crystalline H2O 


13% (40) 


3% (100) 


13% (160) 


12% (150) 


3% (50) 


19% (200) 


Grundv&Schmitt(1998) 


Pure NHi 






2% (120) 








Schmitt etal. (1998) 


NH3:H20(1:99) 




13% (150) 






11% (300) 




T. Roush (Priv. comm) 


Pure CH4 




6% (800) 


5% (800) 


5% (800) 






Ouirico & Schmitt (1997) 


Amorphous C 


71% (10) 


71% (10) 


71% (10) 


71% (10) 


71% (10) 


58% (10) 


Zubkoetal. (1996) 


Titan tholin 


3% (10) 


2% (10) 


2% (10) 


3% (10) 


0% (10) 


1% (10) 


Khare et al. (1984) 


C2H6 








0% (20) 




7% (600) 


Ouirico & Schmitt (1997) 












9% (800) 


13% (800) 


J. Cook (Priv. comm.) 


Reduced 


1.78 


1.29 


1.32 


1.35 


1.64 


1.23 





improve our knowledge about the presence of the various ices 
(such as methane and ammonia) and related chemical products 
and processes on the surface of this object. 



4.7. Surface composition variations over time? 

We saw that the neutral, featureless (within the noise level) vis- 
ible spectral properties are consistent over a timespan of about 
4 years. It is almost impossible to phase the different spectra 
with the rotational ligthcurves available because they were not 
observed close enough in time. However, as Orcus does not 
show dramatic lightcurve effects (see Sect. [U and is expected 
to be a spherical object, it is reasonable to think that the surface 
properties are rather homogeneous over the surface in the op- 
tical. In th e near-infrared, acco rding to the most recent, higher 
SNR data ("Baruc ci et al.ll2008l and this work), crystalline wa- 
ter ice was always present (and in the same roughly propor- 
tions as in the indicative Hapke modeling). A 2.2fim feature was 
also systematically shown to be there in the most recent works, 
with slightly different structures, although we believe they are 
roughly consistent within the noise. Until significantly higher 
SNR data are available and show more subtle changes (as for in- 
stance Charon), we can conclude that Orcus' spectral properties 
are rather homogeneous over its surface in the near-infrared at 
the currently available resolution. 



4.8. Contribution of tine sateiiite 

The known properties of the Orcus binary system are summa- 
rized i n Table [1] Because of the small separation of the pair 
(0.26" INo11 et al.f 2008). it is impossible to separate the contribu- 
tion of both components in our ground-based observations. The 
color measurements performed with the Hubble Spa ce Telescope 

tm which the binary system can be resolved) bv lBrown et alJ 
10) show that Vanth displays significantly redder visible and 
near-infrared colors than Orcus. The contribution of the satellite 
to the total flux is 8% in the visible and a bit more in the near- 
infrared. Thus, we do not expect large variations in the water ice 
content of the main body. Ammonia, methane, and their irradi- 
ation products, if confirmed, should be present at the surface of 
the primary object only, as Vanth is probably too small to have 
retained volatile ices at its surface. Given that the satellite's near- 
infrared color inde x is redder than that of Orcus (as measured by 
iBrown et alj|201(j) . we can expect that the satellite will have a 
lower water ice content, if any. 



5. Comparison with other water-bearing KBOs of 
interest 

In this section we look at the surface expression of water (more 
specifically in its crystalline form) and additional traces of 
volatiles on some other large-intermediate water bearing objects. 
We try to point out common properties and differences, and in 
the case of Quaoar, infer general clues about its possible history 
and interior. The compared properties are summarized in Table 
□ 

5.1. (50000) Quaoar 

The KBOs (90482) Orcus and (50000) Quaoar seem to share the 
same type of near-infrared spectrum: dominated by crystalline 
water ice, with small amounts of volatiles ( Schaller & Brownl 
2007a). Beside s, they are both interm ediate-size KBOs with 
similar albedo dStansberrv et al.l [20081) : their study could thus 
shed light on the possible transition between large methane- 
dominated objects (like Pluto, Eris, Sedna and Makemake) 
and smaller water-ice-dominated ones. Nonetheless, these 
two objects might be very different when it comes to their 
thermal and dynamical history. Moreover, their visible spectral 
properties differ: Quaoar, with its ~28%/100nm slope, is part 
of the objects hosting ultra-red material (>25%/100nm. Oewittl 
2002.) on its surface. Table Q and [8] show a comparison between 
various physical and dynamical properties of Orcus and Quaoar. 

In this table, the maximum surface temperature is computed 
from a thermal balance for the whole surface including different 
energy inputs such as: the incoming solar energy, 

^ (I - Jl)^cos(^), (1) 

(with y[ the Bond albedo, C© the solar constant, dn the he- 
liocentric distance and ^ the local zenith angle), the thermal 
emission scrT'^ (with s the thermal emissivity, cr the Stefan- 
Boltzmann constant and T the temperature), the radial heat 
flux and lateral heat fluxes, considering a thermal conductiv- 
ity of about 10"^ W.m^'.K"'. This thermal balance is com- 
puted for about 50 orbits with a full 3D thermal evolution 
model to get the resulting temperatures (see Appendix B of 
iGuilbert-Lepoutre et al. 1^2010). The insulation per orbit (in J), is 
computed with the formula (Prialnik et aL.2004) 



2VGM0 Vfl(l -e^)' 

with Lq and Mq the solar luminosity and mass respectively, G 
the gravitational constant and R the object's radius. Depending 
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Table 7. Comparison of orbital and physic al pr operties of so me intermediate-large size water-be aring objects of interest. 
References are: (1) Stansberrv et al.' ('2008'), (2) Sicardv et al. ('2006') and references therein, (3) Brown et al. ('2010'), (4) 
[Brown & Trujillo (2004), (5) Brown (2008), (6) Brucker et al. (2009), (7) com puted f rom Buie et al. (2010), (8) Rabinowitz etall 
(2006), (9) Cook et al.i (I2009h . (101 iGuilbert et alJ d2009h . (111 LSchaller & Brownl d2007ah . (121 iDalle Ore et alJ d2009l) . (13) 
[Merlin et ah (2010a) and references therein. 



Parameter 


Orcus 


Charon 


(208996) 2003 AZ84 


Quaoar 


Haumea 


Dynamical Class^ 


3:2 


3:2 


3:2 


Classical 


Classical 


Orbit inclination (°^ 


20.6 


17.1 


13.6 


8.0 


28.2 


Semi-major axis (AU^ 


39.2 


39.6 


39.4 


43.5 


43.0 




30.2 


29.7 


32.3 


41.7 


34.5 


Q(AU:i^ 


48.1 


49.6 


46.5 


45.2 


51.5 


r (AUj^ 


4 / 






A'i 


C 1 

J t 


Diameter (km) 


946±70*" 


1207±3'"* 


686±95<" 


850±200*" 


1150±200*'* 




940±70<3' 






910±100<'" 


~ 2000 X 1500 X 1000<^' 










1260±190<'" 




Pv (%)! 


19.7±3*'* 


~ 40<" 


12.3±4(" 


19.9±10*" 


80±15<'* 




28±4<3' 






17.2±5<« 


60-80"*' 










9.2±4(''> 




Vis. slope (%/100nm:i^ 


~ 2 




~ 3 


~ 28 


~0 


Dominant ice at surface 


H2O (cryst) 


H2O (cryst) 


H2O (cryst) 


H2O (cryst) 


H2O (cryst) 


Volatiles and other 


NHJ+CzHfi 


NH++CH4+NH3:H20<''> 


CH3OH ?<'"> 


CH4+C2H*'" 






NH3:H20+CH4 


NH3:H20,NH3"^' 




CH4+C2H6+N^'^' 






NH3-HCH4 


NH40H<'3) 









° Using classification bv lGladman et alj ( 12008 1) 

* with respect to ecliptic plane, from the Minor Planet Center 
From the Minor Planet Center 

Perihelion from the Minor Planet Center 
Aphelion from the Minor Planet Center 
^ Approximate current heliocentric distance 

* Visible geometric albedo 

'' Slope of the visible spectrum from lFomasier et aP (l2009h and references therein, except for Orcus: this paper 



on the radius value we use for Quaoar, two different values are 
obtained for Q. The values for Orcus and Quaoar (and surface 
temperature) are consistent, although their spectral response to 
solar illumination is very different. Indeed, despite the similar 
visual albedo and near infrared spectral properties, Orcus' 
visible spectrum is neutral, whereas Quaoar's spectrum has a 
red slope; they most probably have different thermal properties, 
in terms of their spectral response to solar illumination, light 
absorption and subsequent heating. 

Quaoar is a classical KBO and as such should have been 
formed in situ, or was moderate ly pushed outward during 
Neptune's migration dGome ^ I2003h . Therefore, it should have 
been formed on a larger timescale than Orcus, implying po- 
tentially less accretional heating. Orcus is a Plutino and might 
have been formed clo ser to the Sun than its current location 
dMalhotralll99lll995h . We can thus suspect that its accretion 
was faster and involved more heat. The early thermal history of 
both objects in this sense might be quite different. Nonetheless, 
because of their rather large sizes and densities, they could 
have both reached temperatures sufficient to melt water ice and 
produce a differentiated internal structure (/. e. after heating by 
short-lived radiogenic elements, the interior should be stratified, 
with a differentiated core and layers enriched/impoverished in 
volatile ices closer from the surface). 

The abundance of radiogenic elements is directly related to 
the dust mass fraction within the body: the denser an icy body 
is, the more heat will be released upon the decay of these ra- 
diogenic elements. Considering Quaoar's high density, we can 



Table 8. Comparison of Orcus and Quaoar additional parame- 
ters. 1: values fr om the Minor Plane t C enter. References are: a : 
lBrownl(l2008b . b: lBrown et"and2010l) . c: lFraser & Brownld2009h 





Parameter 


Orcus 


Quaoar 




Dynamical class 


Plutino 


Classical KBO 


1 


Semi-major axis' (AU) 


39.163 


43.467 





Eccentricity' 
P (yrs) 


0.227 
245.5 


0.040 
286.9 


& 

0- 


Density (g.cm"') 


1.9±0.4'' 
1.5±0.3'' 


>2.8'' 




T,„„-/ max. (K) 


40-50 


42-46 


Phy 


Insulation per orb. (J) 


-2.0x10^' 


1.5x10-' 
3.3x10-' 


CO 


Brightness 




0.6" 


fraction (%) 


8" 



assess that it might have suffered from severe thermal modifica- 
tions due to the decay of long-lived radiogenic elements such as 
'^"K, -^^Th, 2^^U, 2^**U. Therefore, the methane observed at its 
surface could have been placed there in a quite recent past, al- 
though the object should be geologically dead for millions years. 
The subsequent surface evolution should have been dominated 
by space-weathering processes. Orcus is less dense than Quaoar: 
a thermal evolution model is applied to test if long-lived radio- 
genic elements could have heated the object in the last million 
years, so that the effects on the surface would still be observable 
today (see Sect.[6]l. 
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5.2. The Haumea collisional family 



5.4. Plutino (208996) 2003 AZi 



•84 



(136108) Haumea is a 2000 km-class fast-rotating elon- 
gated dwarf planet from the classical Kuiper Belt with the 
highest water ice surface content curre ntly known among 



KBO s and a neutral visi ble spectrum (see iPinilla-Alonso et al 
^009"; 'M erlin et alj 1200 7, and references therein). i Brown et™ 



(2007b) identified several objects with similar spectral prop- 
erties in a dynamical area that is close to Haumea. Because 
Haumea's unusual rotational properties are compatible with 
a giant impact, it was soon suggested that these objects 
are members of the first collisional family identified in 
the Kuiper Belt. Repeated spe ctral a nd photometric studies 
(Snodgrass et al. 
2008;' 



2008; Rabinowitz et al. 



2009; BarkumeetaL 
Pinilla-Al onso et al.i i20()8t ISchaller & Browni i2007bc 



Pinilla-Alonso et al. I l2007h aimed at refining the family proper- 
ties and members list. 

From these studies, we see that despite the presence of water 
ice and the neutral visible spectrum, the other properties of the 
Haumea collisional family differ from those of Orcus. Haumea 
itself has a much higher albedo (0.6-0.84 Rabino witz et alj 
20061: IStansberrv e t al. 2008), and density (2600-3340 kg m'^ 



Rabinowitz et al. 2006) and much deeper water ice absorption 
bands. Pinilla-Alonso et al. (2009) showed that it has a spec- 
trally homogeneous surface best described by a simple intimate 
mixture of amorphous and crystalline water ice and that it is 
probably depleted of any other ices. As a consequence, cryo- 
volcanism cannot be invoked as possible a re-surfacing process, 
whereas it might play a role in Orcus' evolution. In this work, 
we also showed that a simple water ice mixture (Haumea-type 
surface) fails to completely reproduce Orcus' spectral behavior 
in the near-infrared (see model 1). 



5.3. Charon 

Pluto's larger satellite Charon is a ~1200 km size body with a 
water ice dominated surface. Charon and Orcus are of similar 
size and density (1.5-1.9 g.cm"'), on a similar orbit, and share 
common spectral properties. However, their albedos are quite 
different, as Charon was measure d to have a visib le geometric 
albedo of ~0.40 (computed from lBuie et alj|2010l) . Crystalline 
water ice is unambiguously dete cted and leads to a surface 
temperature estimate of 40-50K (ICook et aP l20p7[). which is 
in the same range as that of Orcus (44K from iBarucci et af] 
l2008l) . Another component absorbing at 2.21/zm is preferably 
attributed to a mmonia-type ices from detailed Hapke modeling 
dMerlin et alj 2010a, and references therein), and it seems that 
the surface is variegated (differences between the Pluto and anti- 
Pluto hemisphere s). Recent high SNR K-band spectral studies 
dCook et alj|2009l) suggested for the first time that ammonium 
could be present at the surface of Charon, as it helps to bet- 
ter model the absorption features they detected around 2.21 and 
2.24/vm. Their best-fit Hapke model includes a mixture of water, 
ammonium, methane, and ammonia hydrate. In this work, we 
tried to fit to our data a similar synthetic spectrum: we obtained 
a reduced of 1 .68 (e.g of lower quality than our best-fit op- 
tions), although the general spectral features of Orcus were cor- 
rectly reproduced. This confirms the similarities of Orcus and 
Charon bulk surface properties. Also. lDesch et al.l (12009) identi- 
fied from their ID-spherical thermal model of Charon that cryo- 
volcanism is the most probable mechanism to provide fresh ices 
on the surface. From our own thermal model (see Sect. |6]), we 
draw a similar conclusion for the case of Orcus. 



(208996) 2003 AZ^4 is, like Orcus, a binary Plutino wandering 
in a 30-50 AU heliocentric region. It is presently located at a 
similar heliocentric distance, ~46 AU (wrt ~48 AU for Orcus). 
(208996) has a quasi-neutral featureless spectr um in the opti- 
cal, w ith a spectral slope of 1.6 ± 0.6 %/lOOn m ( Fornasier et al.l 
l2009l and references therein). Guilb ert et al. ( 20091) confirmed 
the presence of water ice (Barkum e et al.l 200 8*) in its crys- 
talline phase from near-infrared spectroscopy, and noted a pos- 
sible absorption feature around 2.3/im. The SNR of the data is 
not sufficient to identify the corresponding absorbing substance. 
However, the general spectral behavior of Plutino (208996) is 
close to that of Orcus, and joint studies deserve to be conducted 
in the future. 



6. A thermal-evolution model of Orcus 

6.1. The model 

We used a thermal-evolution model in order to evaluate the pos- 
sibility for Orcus to undergo some internal activity events such 
as cryovolcanism in a recent past, so that the effects might still 
be observable at the surface. From input orbit, size, albedo, den- 
sity, porosity, formation delay and type of radiogenic elements, 
we simulated the evolution of temperature over the age of the 
Solar System for the whole object. In this model, the heat dif- 
fusion equation (Eq. [3]) is solved using a fully 3D scheme, con- 
sidering internal heating due to the crystallization of amorphous 
water ice, as well as heating due to the decay of several short- 
and long-lived radiogenic elements. We therefore need to solve 
the equation 



dT 



PhuikC -z- + div(-K grad T) = S, 
at 



(3) 



where T is the temperature distribution we need to determine, 
Pbuik the object's bulk density, c the material heat capacity, k its 
effective thermal conductivity and S - Qiad + Qcrvsi the two heat 
sources. The radiogenic heat source is described by 



rad 



y, Qd Xrad Hrad— exp ( — j , 
Trad \Tradl 



(4) 



with Qd the dust bulk density, Xrad the initial mass fraction of 
a particular radiogenic element in the dust, Hrad the heat re- 
leased by this radiogenic element per unit mass upon decay. 
Trad its decay time, and t the time. We consider different short- 
and long-lived radiogenic elements, presented in Table |9] using 
data adapted from Castillo-Rogez et al. (2007, and references 
therein) for ordinary chondrites. 

The heat source due to the crystallization of amorphous wa- 
ter ice is described by 



M T) Qa //„, 



(5) 



with Qa the amorphous water ice bulk density. This phase tran- 
sition releases a latent heat H„c, whose v alue is 9x10"* J .kg"' 
(lKlingeiill98l1) . with a rate measured by ISchmitt et al.l (Il989h 
of: 



A{T) = 1.05X 10'^ e- 



5370/r ^-1 



(6) 



The heat diffusion equation is then expanded on the 
real spherical harmonics basis. The resulting equation is nu- 
merically solved using an implicit stable Crank Nicolson 
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Table 9. Physical parameters of various short-period (SP) and 
long-period (LP) radiogenic elements considered in the model. 
Values adapted from Castillo-Rogez et al. (2007). 





Element 




T * 

(years) 


J.kg-' 




^"Al 


6.1x10-' 


1.05x10" 


4.84x10'^ 


SP 




2.5x10^^/-** 


2.16x10" 


5.04x10'^ 






2.8x10-'* 


5.34x10" 


4.55x10'^ 




40 


1.2x10-" 


1.80x10" 


1.66x10'^ 


LP 


232-pjj 


6.0x10"^ 


2.02x10'° 


1.68x10'^ 


235U 


9.0x10-' 


1.02x10' 


1.82x10'-'' 




238u 


2.9x10-** 


6.45x10' 


1.92x10'-'' 



" Initial mass fraction of the radiogenic element 
* Decay time 

Heat released by this radiogenic element per unit mass upon decay 



scheme. A full description o f this model can be found in 
iGuilbert-Lepoutre et alj(l2010h . 

Several physical parameters have to be initialized before 
simulating the thermal evolution of Orcus. We first assume that 
the object is on its current orbit for the whole simulation and is 
placed there right after its formation time. This does not reflect 
the true dynamical evolution of Orcus, but might not strongly 
affect the resulting general thermal behavior (as detailed in the 
next paragraphs). Orc us is then assumed to have a 475 km ra- 
dius, a 19.9% albedo (Stansberrv et a l.ll2008h . and a bulk density 
of 1900 kg. m'^ ([B rown 2008).' By assuming a residual poros- 
ity of 10% ( Durham et al.l l2005). this density leads to a mixture 
made of 77% of dust and 23% of water ice (mass fractions), 
the dust being homogeneously distributed within the ice matrix. 
The resulting thermal properties are: the thermal conductivity 
K = 0.18 Wm-'K-i, the heat capacity c = lO^Jkg-'K-'. 

Since heating due to the accumulation of gravitational poten- 
tial en ergy is usually not imp ortant, even for large KBOs such as 
Pluto dMcKinnon et al.ll997h . we neglect this heat source in this 
model. We also assume that the objects undergo a cold accretion: 
the accretional heating is thus neglected, although it could be a 
very important heat source in the case of Orcus. 

6.2. Results 

6.2.1 . Effect of short-lived radiogenic elements 

The formation time of KBOs is of critical importance when it 
comes to the study of their differentiation due to heating by 
short-lived radiogenic elements. These elements indeed decay 
during the accretion period, thus providing much less potential 
heat when the bodies are formed (and our simulations start). 
There is still a wide range of available values for this formation 
time, depending on models. For example, Weidenschilling 
( |2004|) suggested that 50 km-radius objects could be formed in 
less than one million years i nside 30 AU. Recent simulations 
performed by iKenvon et al.l (|2008|) show that 1000 km-radius 
bodies could be formed in 5-10x10* yrs within the 20-25 AU 
region. 

With the aim of constraining Orcus' state of differentiation 
(primordial and during its lifetime), we applied the model con- 
sidering various formation delays as input. The initial tempera- 
ture is 30 K (cold accretion), and the radiogenic elements initial 
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Fig. 5. Evolution of Orcus' central temperature (black circles) as 
a function of the simulated time in the case of a formation delay 
of 1 million years and a high mass fraction of radiogenic element 
''"Fe. Triangles and square curves correspond to the temperature 
at a depth of 1km and 500m resp. Top graph zooms on the first 
~10 Myr of simulation. 



mass fraction is calculated using the following formula, to ac- 
count for their decay during accretion 

X,-ad(tF) = X,ad{Q) e-"l'"-\ (7) 

with tp the considered formation delay, and Xradi^) found in 
Table |9l All SP and LP elements listed there are taken as in- 
put, and we distinguish the case of the high and low ^''Fe mass 
fraction. Our simulations over the age of the Solar System 
(4.5x10'^ yrs) show that Orcus could have a physically differ- 
entiated structure -including a melted core, crystalline water ice 
mantle, and a pristine surface layer composed of amorphous ice 
and dust- if formed in less than about 5x10* yrs. The melting 
point of water is indeed achieved within the body in the first mil- 
lions years of the simulations. If we assume that ammonia is also 
present in the mixture, then this formation time limit is slightly 
increased to about 7x10* years, since the presence of ammonia 
depresses the melting poi nt from a few K to almost 100 K de- 
pending on its abundance (lKarge]|ll992L and references therein). 

Figure |5] shows the temperature evolution inside Orcus at 
various depths (center, Ikm-deep, 500m-deep) considering the 
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Fig. 6. Temperature (K) as a function of the object's depth for 
various evolution times. This plot illustrates how the outside 
medium cools the surface layers deeper with time (for SP and 
LP radiogenic elements, high *"Fe mass fraction.) 
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Fig. 7. Central temperature (K) as a function of the simulated 
time, showing the different contributions from SP+LP radio- 
genic elements, and LP only. 



highest mass fraction of short-lived Fe and a formation delay 
of 1 Myr Although this formation delay might seem too short 
for such a large object, the corresponding results show the great- 
est temperature variations and the steepest thermal gradients ex- 
pected, to put an upper limit to all thermal processes. The radio- 
genic heating is homogeneous inside the body, but the surface is 
in thermal equilibrium with the external medium. This results in 
the propagation of a cold wave from the surface toward deeper 
layers. As a consequence, the curves corresponding to the 1km- 
depth and 500m-depth show a maximum temperature after a few 
millions years of the simulations: the body then starts to cool 
slowly. 

Another consequence is that a layer of amorphous water ice 
is maintained at the surface of the object: the Ikm-deep layer 
temperature reaches the water crystallization threshold at some 
point, while the 500m-deep layer does not complete the crys- 
tallization process. The propagation of the cold wave from the 
surface to the interior is illustrated by Fig. |6l which shows the 
radial temperature profile across the body for different simulated 
times. The thickness of this pristine material layer depends on 
the formation delay that we consider: the shorter the delay, the 
thinner the layer. We find that the thickness could be as low as 
several hundreds of meters (~500m for the 1 Myr delay) or as 
thick as a few tens of kilometers (~10km for the 10 Myr delay). 
Most of the volume is therefore made of crystalline water ice. 
Our simulations also show that the threshold for the solid-liquid 
phase transition of water is reached within most of the volume (if 
formed fast enough): the high temperatures needed to keep liq- 
uid water inside Orcus are moreover maintained for long periods 
such as 2x10^ years. 

6.2.2. Effect of long-lived radiogenic elements 

What would happen if Orcus was formed in more than 10 Myr? 
In this case, the effect of short-lived radiogenic elements would 
be nonexistent. However, radiogenic heating due to the decay of 
40k, 232Th, 235u, and ^^^u would still be inevitable. This would 



even be of the greatest importance for Orcus, whose thermal 
timescale. 



exceeds the age of the Solar System (assuming the previous ther- 
mal physical parameters). Nonetheless, the value of the thermal 
conductivity is critical for this calculation, and a variation of one 
magnitude could reduce the cooling time to less than 10^ years. 

We therefore applied our model including only the long-lived 
radiogenic elements "'"K, ^^^Th, '^^^\J, and '^^^\J. The evolution 
of the body's central temperature is shown on Fig.|7] compared 
to the evolution under the effect of both long- and short-lived 
elements. Even without the high heating power of short-lived 
elements, the internal temperature increases, reaching the melt- 
ing point of water ice within the first billion years. Please keep 
in mind that our mixture contains a high fraction of dust that 
might not reflect the true composition of Orcus. In this case 
Orcus suffers from considerable heating, while it could actually 
be much more moderate. Our simulations show that in the case 
of an important heating by long-lived radiogenic elements, the 
melting point of water ice could be reached within about 60% of 
the body, while the melting point of an ammonia-water mixture 
would be reached within 80% of the body. This is illustrated in 
Fig. [8] which shows radial temperature profiles for different sim- 
ulated times. As in the cases considered in the previous section, 
a layer of amorphous water ice is maintained at the surface of 
Orcus. Its thickness is about 5-6 km. 

6.3. Implications 
6.3.1. Limitations 

Our model has several limitations that need to be kept in 
mind when interpreting the results. The main issue is that 
the liquid phase or differentiation is not accounted for: the 
model is not suitable anymore when the temperature threshold 
for solid-liquid phase transition is reached. The results we 
provide should therefore be regarded as qualitative evolutionary 
trends rather than a quantitative description of Orcus' thermal 
evolution. Accounting for sources of substantial cooling such 
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Fig. 8. Temperature (K) as a function of the object's depth for 
various evolution times (the lower limit is the case of radio- 
genic heat source from long-lived elements only). Horizontal 
lines show the temperature above which H2O or H20-HNH3 are 
in the liquid phase. 



as convection movements would probably induce moderate 
or even low temperature at the present time. Moreover, we 
considered a water ice and dust mixture which contains a high 
fraction of dust, which in turn induces a high heating power by 
radiogenic elements. This dust fraction might actually be lower, 
and some thermal properties such as the thermal conductivity 
might also be reduced by a large factor in reality. However, even 
when considering low heating, our simulations show that Orcus 
should have experienced a cryovolcanism episode in its life, 
especially if ammonia is present in its bulk composition. 

We did not take into account the possible contribution from 
the satellite in our simulations, as it is currently not possible 
to constrain it. However, if the binary system was created by 
a giant impact, we can expect that a large fraction of the main 
body's surface layers was removed, thus exposing interior ma- 
terial. Furthermore, the heat from the impact should contribute 
to a temporal increasing of Orcus' upper layers temperature. 
Tidal heating was not accounted for either, and might play a role 
in Orcus thermal evolution if the sat ellite is found to b e large 
enough (current rough estimates from lBrown et alj|2010l give a 
size of the satellite of 1/3 to half the primary). 

6.3.2. A past cryovolcanic event 

Under the effect of an efficient heating due to radiogenic ele- 
ments, a liquid phase (water -i-ammonia, if ammonia is present) 
is produced, and could have reached layers close to the surface, 
eventually freezing, and/or even propagating thro ugh cracks 
withi n the ice matrix, thus erupting at the surface (IStevensonI 
11982 "). This process could be even more efficient with methane, 
which can trigger explosive cryovolcanism (Kargel 1992, 1995). 
Our spectral data show a compatibility with traces of ammo- 
nia and methane at the surface of Orcus, which means that the 
interior of the object most probably contains ammonia. There 
is nonetheless no easy way to infer the internal abundance of 
ammonia from the surface composition, since the latter would 
not reflect the internal composition, due to the chemical modi- 
fications induced by space weathering processes. As described 



by iDesch et"an (l2009h . the presence of only a few percent of 
ammonia in the interior (by weight, relative to water) can effi- 
ciently lower the viscosity of the ammonia-water mixture once 
is has melted. Please note however that large amounts of ammo- 
nia might be incompat ible with a high du st/ice ratio such as the 
one used on our model dWong et al.l2008h . Our simulations show 
that we cannot rule out a cryovolcanic episode in Orcus' past 
life, which probably supplied most of the observed crystalline 
water ice of the surface. Surface erosion with impact cratering 
could also contribute to bring to the surface some fresh material, 
including crystalline water ice, ammonia, and methane, refrozen 
in su bsurface layers, but to a much lesser extent (as computed by 
ICook et al...2007. for Charon). 

6.3.3. Origin and surviving of crystalline water ice 

The initial state of water ice in the Solar System is still a matter 
of debate: 

- If water ice was initially crystalline when Orcus formed, then 
it could have stayed so over the age of the Solar System in 
the object's interior Space weathering processes would have 
induced its progressive amorphization at the surface (partial, 
as we will see below, or complete). 

- If water ice was initially amorphous when Orcus formed (as 
we assumed), then our simulations show that it would have 
been rapidly crystallized compared to the age of the Solar 
System for most of the objects' volume. A layer of pristine 
amorphous ice (a few kilometers thick in the canonical case) 
would have been maintained at the surface. In this case, the 
crystalline water ice would only appear at the surface if a 
supply mechanism was at work at some point of the his- 
tory, or if the erosion by impacts removed the amorphous 
ice over few kilometers (which is currently believed to be 
an inefficient process to explain the detectable quantities of 
crystalline water ice on Orcu s, other KBOs an d giant plan- 
ets' satellites, as computed bv lCook et alj|2007l for Charon). 
However, a more detailed modeling of the impact environ- 
ment of Orcus needs to be conducted before completely rul- 
ing out the resurfacing effects of impacts. 

In our particular case, the detection of large amounts of 
crystalline water ice at the surface of Orcus could indicate that 
heating from radiogenic decay and the crystallization of amor- 
phous ice had been an efficient process, thus bringing the crys- 
tallization front close to the surface, as shown by our thermal 
evolution model. But most importantly, our results show that 
Orcus' thermal history is compatible with a past cryovolcanic 
event, which is an efficient supply mechanism of fresh ices to 
the surface. We can therefore envision that most of the crys- 
talline water ice we currently observe are the actual leftovers 
from suc h a geological even t. Indeed, recent laboratory experi- 
ments by IZheng et aP (1200^ ). show that the 1 .bSfim crystalline 
water ice absorption band is able to survive ion irradiation doses 
up to 160+30 eV for T>30K. This phenomenon is temperature- 
dependent, and becomes even more efficient at higher temper- 
atures. This means that the amorphization of the crystalline ice 
is not complete, and starting at 50K, there is a balance dom- 
inated by ther mal recrystall i zation . From the irradiation doses 
determined by ICooper et al.l (l2003b for the optical thickness of 
0.1 mm {i.e. the layers probed by near infrared spectroscopy), 
Zheng et al. computed for a KBO, that the 1.65/vm band should 
be able to survive during about 6 billion yrs, e.g. longer than the 
Solar System age. If this result is confirmed, it means that the 
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40-50K surface of Orcus still holds traces of a past cryovolcanic 
event in the form of detectable crystalline water ice, and that ad- 
ditional resurfacing processes are not necessarily needed. 

7. Summary 

We provided new 0.4-2.35 /ym reflectance spectra and near- 
infrared photometry of Plutino (90482) Orcus. Our visible 
spectrum is featureless (within the SNR of the data) with a 
neutral slope (~2%/100nm). These results are consistent with 
previously published works and point (within the noise level) to 
a homogeneous surface at least in the optical. Our near-infrared 
spectra confirm the presence of water ice in its crystalline state 
(thanks to the 1.65/im absorption band) at the surface of Orcus, 
as well as the presence of a 2.2yum feature, which cannot be 
readily identified owing to both the detection level and the 
multiple candidates. This band might show a double structure, 
which deserves to be further investigated from significantly 
higher SNR studies than currently available. These general 
spectral properties are shared with Pluto's satellite Charon and 
Plutino (208996) 2003 KLu- 

Radiative transfer Hapke models show that a simple mixture 
of water ice (amorphous and crystalline) is not sufficient to 
properly describe the data in the H and K band. Orcus therefore 
does not have a Haumea-type surface, although H2O should 
be the dominant ice. Additional absorbents should be invoked, 
which may be volatiles and their possible irradiation products. 
Further Hapke modeling shows that the 2.2yL(m region and 
overall spectrum from 0.4 to 2.35 [im are best described by 
an intimate mixture of ammonium (NHp and traces of ethane 
(CaHg), which are most probably solar irradiation products, 
in addition to water ice (mostly in its crystalline state). The 
presence of ammonium even helps to better reproduce the 
visible spectrum, an effect which was not expected a priori. 
We find also that our data are compatible with Hapke models 
of an equivalent (slightly lower) statistical level including a 
mixture of water ice, methane, and ammonia (pure of diluted); 
the latter are the parent species of ethane and ammonium. 
Methane is particularly interesting to reproduce the blue slope 
longward ~2.3;um, which might also be the shortward wing 
of an additional absorption feature. Focused studies should be 
conducted in close connection with laboratory investigations to 
better understand the methane and ammonia chemistry. 

In this work, we showed that the Orcus spectrum could 
be compatible with the pre sence of volatile ice s on its surface 
(as also hypothesized by Barucci et al.l I2008), and their in- 
volatile irradiation products (first invoked here). Our findings 
merits confirmation by higher SNR data, and we believe that 
Plutino (208996) 2003 AZ84 and Orcus deserve joint detailed 
studies, to be also connected with Charon properties , as they 
share common spectral properties. ISchaller & BrowrJ (l2007bh 
developed a simple model of atmospheric escape to predict how 
a given KBO can preserve a primordial surface inventory of 
CO, N2 and CH4. Orcus, Charon and Quaoar are located in a 
relative narrow region in their Fig. 1, showing the equivalent 
temperature (K) as a function of the object's diameter (km). 
Orcus is located in the "all volatiles lost" region, where methane 
and nitrogen ices in particular should have completely escaped 
from the surface. That our Orcus data could be compatible 
with the presence of volatiles (here CH4) would imply that (i) 
such an atmospheric model is too simple to describe the actual 
volatile content of a particular object (as it was emphasized 



by the authors in their conclusions); (ii) Orcus might be larger 
than the size we currently know from far-IR data (and Herschel 
Space Observatory observations might allow us to answer this 
question, although from the Schaller and Brown model, a factor 
of ~1.7 would be needed); (iii) there is a supply of volatile 
ices from the interior of the object; (iv) the actual situation is a 
combination of the previous points. 

We investigated the thermal properties of Orcus over its 
lifetime to understand how crystalline water ice can be produced 
and test if geological events such as cryovolcanism can be 
possible and provide fresh ices to the surface. We used a 3D 
thermal evolution model by iGuilbert-Lepoutre et al.l (12010) that 
simulates the temperature of the whole object as a function of 
its size, density, porosity, orbit, formation delay and radiogenic 
elements primordial inventory (^''Al, ''°Fe, ^^Mn, "^''K, ^''^Th, 
235jj^ -pjjg ^^jjj- assumed composition is dust (77%) 

plus amorphous water ice (23%), the relative mass fractions 
being derived from the input 10% porosity. Various initial 
configurations have been simulated with different formation 
delays (from 1 to over 10 Myr) including different amounts of 
short-lived radiogenic elements (with the case no short-lived 
elements at all) to evaluate the differentiation state of Orcus. Our 
simulations show that a layer of amorphous water ice of a few 
kilometers thick would be maintained at the surface of the body, 
while the remaining of the volume is occupied by crystalline 
water ice. For all cases we considered, the heating induced by 
radiogenic elements leads to the melting of the core. This effect 
is strengthened if ammonia is present inside the body, because 
it reduces the melting point temperature of the mixture and 
allows for the production of a liquid phase even under moderate 
heating. This liquid phase could easily reach layers close to the 
surface through cracks in the ice matrix (if not directly produced 
in subsurface layers in the more extreme cases), and therefore 
make Orcus and intermediate to large water-bearing KBOs the 
subject of potential astrobiological interest. 

As a consequence, even if it is premature to suggest that 
Orcus did undergo a recent event of cryovolcanism, our results 
show that this event did take place at some point of Orcus' his- 
tory and most probably supplied cryst alline water ic e to the sur- 
face. Recent laboratory experiments bv lZheng et al.l (12009) show 
that the ion-irradiation-induced amorphization of crystalline wa- 
ter ice on a KBO at ~40K can be incomplete and be efficiently 
balanced by re-crystallisation (starting at 50K), so that crys- 
talline water ice can survive and still be detectable at the surface 
over timescales longer than the Solar System age. Both from the 
latter results and our model, we can envision that some of the 
crystalline ice we currently observe at the surface of Orcus is 
actually the remnant of a past cryovolcanic event. In addition, 
in the most probable case of our simulations, a large volume of 
volatiles are currently present in Orcus subsurface layer, await- 
ing to be revealed by future surface alterations. 
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